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An rf superconducting quantum interference device (SQUID) array in an alternating magnetic 
field is investigated with respect to its effective magnetic permeability, within the effective medium 
approximation. This system acts as an inherently nonlinear magnetic metamaterial, leading to 
negative magnetic response, and thus negative permeability, above the resonance frequency of the 
individual SQUIDs. Moreover, the permeability exhibits oscillatory behavior at low field intensities, 
allowing its tuning by a slight change of the intensity of the applied field. 
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An rf superconducting quantum interference device 
(SQUID) consists of a superconducting ring interrupted 
by a Josephson junction (JJ)^ When driven by an alter- 
nating magnetic field, the induced supercurrents around 
the ring are determined by the JJ through the celebrated 
Josephson relations. This system exhibits rich nonlinear 
behavior, including chaotic effects^ Recently, quantum rf 
SQUIDs have attracted great attention, since they consti- 
tute essential elements for quantum computing^ In this 
direction, rf SQUIDs with one or more zero and/or 7r fer- 
romagnetic JJs have been constructed. 4 In this Letter we 
show that rf SQUIDs may serve as constitual elements for 
nonlinear magnetic metamaterials (MMs), i.e., artificial, 
composite, inherently non-magnetic media with (positive 
or negative) magnetic response at microwave frequencies. 

Classical MMs are routinely fabricated with regular ar- 
rays of split-ring resonators (SRRs), with operating fre- 
quencies up to the optical ranged Moreover, MMs with 
negative magnetic response can be combined with plas- 
monic wires that exhibit negative permittivity, produc- 
ing thus left-handed (LH) metamaterials characterized 
by negative refraction index. Superconducting SRRs 
promise severe reduction of losses, which constrain the 
evanescent wave amplification in these materials^ Thus, 
metamaterials involving superconducting SRRs and/or 
wires have been recently demonstrated experimentally. 7 
The effect of incorporating superconductors in LH trans- 
mission lines has been also studied, 8 Naturally, the theory 
of metamaterials has been extended to account for non- 
linear effects^ ' 10 ! 11 ' 12 ! 13 Nonlinear MMs support several 
types of interesting excitations, e.g., magnetic domain 
wallsjii discrete breathers^ and evnelope solitonsi^ 3 - 
Regular arrays of rf SQUIDs offer an alternative for the 
construction of nonlinear MMs due to the nonlinearity of 
the J J. 

Very much like the SRR, the rf SQUID (Fig. 1(b)) is 
a resonant nonlinear oscillator, and similarly it responds 
in a manner analogous to a magnetic "atom" in a time- 




FIG. I: Schematic drawing of the SQUID array, along with 
the equivalent circuit for an rf SQUID in external flux & e xt- 



varying magnetic field with appropriate polarization, ex- 
hibiting a resonant magnetic response at a particular fre- 
quency. The SRRs are equivalently RLC circuits in sc- 
ries, featuring a resistance R, a capacitance C and an in- 
ductance L, working as small dipoles. In turn, adopting 
the resistively and capacitively shunted junction (RCSJ) 
model for the JJji the rf SQUIDs are not dipoles but, 
instead, they feature an inductance L in series with an 
ideal Josephson element (i.e., for which I = I c sin <f>, with 
(j> the Josephson phase), shunted by a capacitor C and 
a resistor R (Fig. 1(c)). However, the fields they pro- 
duce are approximatelly those of small dipoles, although 
quantitatively they are affected by flux quantization in 
superconducting loops. Consider an rf SQUID with loop 
area S = ira 2 (radius a), in a magnetic field of amplitude 
H e o, frequency to, and intensity H ext — H e o cos(ujt) per- 
pendicular to its plane (t is the time variable). The field 
generates a flux <3> ea;t = <I> e o cos(wi) threading the SQUID 
loop, with $eo = fioSHeo, and fig the permeability of the 
vacuum. The flux $ trapped in the SQUID ring is given 
(in normalized variables) by 



/ = u 



pi, 



(1) 



where / = $/$ , fext = ^ext/^o, i = I/h, P = Pl/2tt = 
LI c /&o, I is the current circulating in the ring, I c is 
the critical current of the J J, L is the inductance of the 




FIG. 2: (Color online) (a) Coefficient D vs. the applied flux 
amplitude f e o, for j3 — 0.15 (red-solid curve); (3 = 0.10 (green- 
dashed curve); /3 — 0.05 (blue-dotted curve), (b) The ampli- 
tude of the flux /o vs. /eo, for fi = 0.9 (red-solid curve), 
fi = 1.1 (green-dashed curve), and 7 = 0.001, (3 = 0.15. The 
black circles and blue diamonds correspond to the numerically 
obtained fo for = 0.9 and 1.1, respectively. 



SQUID ring, and $0 is the flux quantum. The dynamics 
of the normalized flux / is governed by the equation 

+ 7 |: +/? sin (2^f) + / = / e;ct , (2) 

where C and R is the capacitance and resistance, respec- 
tively, of the J J, 7 = Lloq/R, t — ujQt, loq — \/LC, and 



fext = / e oCOs(fiT), 



(3) 



with / e o = $ e o/ < l , Oj and f2 = oj/u)o. The small parameter 
7 actually represents all of the dissipation coupled to the 
rf SQUID. 

An approximate solution for Eq. ([2]) may be obtained 
for ft close to the SQUID resonance frequency (ft ~ 1) 
in the non-hysteretic regime (3l < 1- Following RefJ^ 
we expand the nonlinear term in Eq. |(3J| in a Fourier - 
Bessel series of the form 

00 1 _ ^ \ n 

/3sin (2tt/) = - V - — -J n (np L ) sin(2Trnf ext ), (4) 

* * 777T 



where J n is the Bessel function of the first kind, of or- 
der n. By substiting Eq. ([3]) in Eq. ([4]) and carrying 
out the Fourier - Bessel expansion of the sine term, one 
needs to retain only the fundamental ft component in the 
expansions— This leads to the simplified expression 



/3sin (27r/)~ J D(/e )cos(fir), 



(5) 



where D(f e0 ) = -2££ 
substitution of Eq. 



i-iy 



= i —^rJn(nf3 L )Ji(2Tmf e0 ). By 
© in Eq. ©, the latter can be 
solved for the flux f — fo cos(ftr + 0) in the loop, with 



fo 



feO-D 



v^ft 2 + (1 - ft 2 ) 2 



6 = tan -1 



-7ft 



1 -ft 2 



,(6) 



where is the phase difference between / and fext- The 
dependence of D and /o on f e o for low field intensity 
is illustrated in Figs. 1(a) and 1(b), respectively. For 
larger f e0 the coefficient D approaches zero still oscillat- 
ing, while /o approaches a straight line with slope de- 
pending on ft and 7. For 7 <C 1 and not very close to the 
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FIG. 3: (Color online) Time-dependence of the flux / (green- 
dashed curves), the applied flux f ex t (blue-dotted curves), and 
the response /3i (red-solid curves), for (5 — 0.15, 7 = 0.001, 
/ e0 = 1.14, and (a) Q, = 0.63; (b) Q = 0.9; (c) fi = 1.1; (d) 
Q = 9.0. 



resonance, 9 ~ 0. It is instructive to express the 7 = 
solution as: 

/ = ±|/o|cos(fir), |/o| = (/ e0 -£>)/|l-0 2 |. (7) 

The plus (minus) sign, corresponding to a phase-shift 
of (tt) of / with respect to f ex t, is obtained for ft < 1 
(ft > 1). Thus, the flux / may be either in-phase (+ 
sign) or in anti-phase (- sign) with f ex t, depending on ft. 
This is confirmed by numerical integration of Eq. 
as shown in Fig. 2, where we plot separately the three 
terms of Eq. ([T]) in time. The quantities /, fext-, and 
fti are shown for two periods T — 2tt/H in each case, 
after they have reached a steady state. For ft < 1 (Figs. 
2(a) and 2(b)), the flux / (green-dashed curves) is in- 
phase with fext (blue-dotted curves), while for ft > 1 
(Figs. 2(c) and 2(d)) the flux / is in anti-phase with 
fext- The other curves (red-solid curves) correspond to 
Pi, the response of the SQUID to the applied flux. Away 
from the resonance, the response is (in absolute value) 
less than (Fig. 2(a), for ft = 0.63) or nearly equal (Fig. 
2(d), for ft = 8.98) to the magnitude of fext- However, 
close to resonance, the response f3i is much larger than 
fext, leading to a much higher flux / (Figs. 2(b) and 
2(c) for ft = 0.9 and ft = 1.1, respectively). Moreover, 
in Fig. 2(c), / is in anti-phase with f ex t, showing thus 
extreme diamagnetic (negative) response. The numeri- 
cally obtained amplitudes fo (depicted as black circles 
for ft = 0.9 and blue diamonds for ft = 1.1 in Fig. 2(b)) 
are in fair agreement with the analytical expression, Eq. 
©. The agreement becomes better for larger / e0 . 

We now consider a planar rf SQUID array consisting of 
identical units (Fig. 1(a)), and forming a lattice of unit- 
cell-side d; the system is placed in a magnetic field H ex t = 
H perpendicular to SQUID plane. If the wavelength of 
H is much larger than d, the array can be treated as an 
effectivelly continuous and homogeneous medium. Then, 
the magnetic induction B in the array plane is 



B = n Q {H + M) = fioHrH, 



(8) 



3 



2.5 



1.5 



■ 



Imm 




0.7ir 



2 4 6 8 10 2 4 6 8 10 



e0 



eO 



FIG. 4: (Color online) Relative permeability fi r vs. / e o, for 
F — 0.01 (red-solid curves), F — 0.02 (green-dashed curves), 
F = 0.03 (blue-dotted curves), and Q, = 0.99 (a); 1.01 (b). 



where M = S I/d 3 is the magnetization induced by the 
current / circulating a SQUID loop, and \i r the relative 
permeability of the array. Introducing M into Eq. ([5]), 
and using Eqs. (P), ©, and (0, we get 



» r = l + F(±\f \/f e0 -l). 



(9) 



where F = Tr 2 (fj,oa/ L)(a/d) 3 . The coefficient F has to 
be very small (F <C 1), so that magnetic interactions be- 
tween individual SQUIDs can be neglected in a first ap- 
proximation. Recall that the plus sign in front of |/o|// e o 
should be taken for 51 < 1, while the minus sign should be 
taken for fi > 1. In Fig. 3 we plot /j, r both for il < 1 (Fig. 
3(a)) and > 1 (Fig. 3(b)), for three different values of 
F. In real arrays, that coefficient could be engineered to 
attain the desired value. In both Figs. 3(a) and 3(b), 
the relative permeability /i r oscillates for low intensity 
fields (low fext), while it tends to a constant at larger 
f ex t. In Fig. 3(a) (Q < 1), the relative permeability fi r 
is always positive, while it increases with increasing F. 
In Fig. 3(b), however, \x r may assume both positive and 
negative values, depending on the value of F. With ap- 
propriate choise of F, it becomes oscillatory around zero 



(green-dashed curve in Fig. 3(b)) allowing tuning from 
positive to negative fi r with a slight change of f ex t- 

In conclusion, we have shown that a planar rf SQUID 
array exhibits large magnetic response close to resonance, 
which may be negative above the resonance frequency, 
leading to effectivelly negative fi r . For low field inten- 
sities (low fext), l^r exhibits oscillatory behavior which 
gradually dissappears for higher f ext . This behavior may 
be exploited to construct a flux-controlled metamaterial 
(as opposed to voltage-controlled metamaterial demon- 
strated in Rcf. 16 ). The physical parameters required for 
the rf SQUIDs giving the dimensionless parameters used 
above are not especially formidable. An rf SQUID with 
L ~ 105 pH, C ~ 80 fF, and I c ~ 3 fiA, would give 
(3 ~ 0.15 (Pl — 0.94). For these parameters, a value of 
the resistance R ~ 3.6 K£l is required in order to have 
7 ~ 10~ 3 , used in the numerical integration of Eq. @. 
However, our results are qualitatively valid for 7 even 
an order of magnitude larger, in which case R ~ 360 tt. 
We note that u>q — u> p /\/]3l, where oj p is the plasma fre- 
quency of the JJ. For the parameters considered above, 
where fir, is slightly less than unity, the frequencies u>q 



and ujp are of the same order. 



However, u> p does not seem 



to have any special role in the microwave response of the 
rf SQUID. Du et al. have studied the quantum version 
of a SQUID array as a LH metamaterial, concluded that 
negative refractive index with low loss may be obtained 
in the quantum regime Consequently, fi r can be neg- 
ative at some specific frequency range. However, their 
corresponding expression for [i r is linear, i.e., it does not 
depend on the amplitude of the applied flux, and thus it 
does not allow flux-tuning. Moreover, experiments with 
SQUID arrays in the quantum regime, where individual 
SQUIDs can be described as two-level systems, are much 
more difficult to realize. 
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